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The magnetic ordering of the hexagonal multiferroic compound YbMnOs has been studied be- 
tween 100 K and 1.5 K by combining neutron powder diffraction, 170 Yb Mossbauer spectroscopy 
and magnetization measurements. The Yb moments of the two crystallographic sites order at two 
different temperatures, the 46 site together with the Mn moments (at Tn —85 K) and the 2a site well 
below (at 3.5 K). The temperature dependences of the Mn and Yb moments are explained within 
a molecular field model, showing that the 4b and 2a sites order via Yb-Mn and Yb-Yb interactions 
respectively. A simple picture taking into account the local Mn environment of the Rare earth R 
(4b) ion is proposed to couple R and Mn orders in hexagonal RMnOs manganites. The nature and 
symmetry of the R-Mn interactions yielding the R order are discussed. 

PACS numbers: 75.25.+z,76.80.+y,75.80.+q,75.30.Kz 



I. INTRODUCTION 

RMnOs manganites, where the R 3+ ions are either Y 
or a rare earth, belong to the multiferroic family, showing 
a coupling between magnetic and electric order param- 
eters. This magneto-electric (M-E) coupling opens the 
possibility of tuning the magnetization by an electric field 
or vice- versa, with potential application for building new 
kinds of electronic devices^. Numerous experiments in 
the RMnOs series probe the key role of their non collinear 
magnetic structures, induced by frustrated magnetic in- 
teractions, in driving such M-E coupling 3 . The M-E cou- 
pling is also connected with a spin-lattice coupling, re- 
cently observed both in structural 4 and dynamical^^ 
properties. Understanding the microscopic characteris- 
tics of the magnetic order and the origin of the magnetic 
interactions is a key issue for tuning their properties. 

The RMn03 compounds are divided in two classes 
according to the ionic radius of the rare earth. Or- 
thorhombic structures are stabilized for large ionic radius 
of the R 3+ ion (R=La-Dy), whereas hexagonal struc- 
tures are stable at small ionic radius (R=Y, Ho-Lu). 
In orthorhombic compounds where magnetic frustration 
arises from competing exchange interactions, ferroelec- 
tric and magnetic orders appear at the same tempera- 
ture (~60 K), which suggests a strong M-E coupling—. In 
hexagonal compounds, magnetic frustration arises from 
the lattice geometry, since the triangular Mn lattice is 
frustrated for antiferromagnetic (AF) first neighbor in- 
teractions. Here, the ferroelectric order occurs around 
900 K, well above the magnetic ordering temperature 
(^80 K). The M-E coupling is evidenced for instance by 
anomalies in the dielectric constant 9 or specific heat^ at 
the magnetic transitions. 

In spite of numerous experiments, the complex mag- 
netic structures of the RMnOs are still not fully under- 
stood. This is especially true for the hexagonal com- 
pounds, where homometric configurations of the Mn mo- 



ments in the triangular Mn lattice are extremely diffi- 
cult to distinguish, and where the R magnetic order is 
intricate, owing to the presence of two different crystal- 
lographic rare earth sites 10 . 

Here we present a detailed study of the magnetic order 
in YbMn03, using neutron diffraction, Mossbauer spec- 
troscopy on the isotope 170 Yb and magnetization mea- 
surements. The magnetic structure of this compound 
has not been studied up to now, although precise predic- 
tions could be made from optical spectroscopy data 11 . 
By combining microscopic and macroscopic probes in a 
large temperature range (1.5 K<T<100 K), we provide 
a complete determination of the magnetic structure. 

The T-dependence of the Mn and Yb moments de- 
rived from our experiments can be explained quantita- 
tively within the molecular field model. This allows us 
to clarify the respective roles of Mn-Mn, R-Mn and R- 
R interactions in this compound. Using these results, 
we discuss the symmetry and the possible origins of the 
magnetic R-Mn interactions, and we propose a simple 
picture, which may hold for the whole hexagonal family. 



II. EXPERIMENTAL DETAILS 

A polycrystalline sample of YbMnOs was prepared 
as described in Refll2l and characterized using X-Ray 
diffraction at room temperature, showing that the sam- 
ple is single phased. Magnetic measurements were per- 
formed with a superconducting quantum interference de- 
vice (SQUID) magnetometer between 1.5 K and 300 K. 

Mossbauer spectroscopy on 170 Yb realises 7-ray reso- 
nant absorption between the ground nuclear state, with 
spin I g =0, and the first excited state, with spin I e =2, 
at an energy of 84.3 keV. We used a neutron activated 
Tm*Bi2 7-ray source, mounted on a triangular veloc- 
ity electromagnetic drive. 170 Yb Mossbauer absorption 
spectra were recorded in the temperature range 1.35 - 
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TABLE I: Atom positions and discrepancy factors at room 
temperature 



50 K. 

The crystal structure was determined by measuring a 
neutron powder diffraction (NPD) pattern at 300 K on 
the high resolution powder diffractometer 3T2 of the Lab- 
oratoire Leon Brillouin (LLB) at Orphee reactor, with an 
incident neutron wavelength A = 1.2253 A. The magnetic 
structure was studied by collecting NPD patterns at sev- 
eral temperatures, between 100 K (above the magnetic 
transition) and 1.5 K, on the neutron diffractometer G4- 
1 of the LLB, with an incident neutron wavelength of 
2.4226 A. Both crystal and magnetic structures were re- 
fined using the Fullprof suited. 



III. CRYSTAL STRUCTURE 

The refined NPD pattern at 300 K is shown in FigtH 
All Bragg reflexions of the pattern can be indexed within 
the hexagonal space group P6scm. The lattice constants 
at 300 K are a = 6.0629(1) A and c = 11.3529(1) A. 

The refined atomic positions are reported in Table H 
They are in good agreement with previous determina- 
tions from X-ray diffractio n 14 ! 15 . As already stated in 
Refll6l the lattice unit cell consists of 6 formula units. 
Each Mn atom is surrounded by oxygen ions forming a 
MnOs bipyramidal structure, with 3 O (2 O4 and one 
O3) ions close to the Mn plane, and two O (Oi and O2) 
ions at the apexes. Corner sharing MnOs bipyramids 
form layers separated along the c-axis by Yb layers. In 
the space group P63C171, Yb occupies the two crystallo- 
graphically inequivalent sites 2a and 46, with trigonal lo- 
cal symmetry around the hexagonal c-axis. The Yb(4b) 
site consists of 4 equivalent atomic positions inside the 
hexagonal unit cell, and the Yb(2a) site, of two equiva- 
lent atomic positions along the c-axis-edges of the cell. 



IV. 



MAGNETIC SUSCEPTIBILITY AND 170 YB 
MOSSBAUER DATA 



The thermal variation of the magnetic susceptibility 
x(T), measured in a field of 20 G, is shown in Fig|2j 
The inverse susceptibility follows a Curie- Weiss law in 
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FIG. 1: (Color online) Observed and Fullprof calculated NPD 
pattern at room temperature. The Bragg reflections (tics), 
and the difference between the observed and calculated pat- 
terns are plotted at the bottom. 
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FIG. 2: (Color online) Inverse magnetic susceptibility in 
YbMnOs- The dashed line is a Curie- Weiss law with 
Me//=6.1/ie and P = — 220 K. Insert: low temperature sus- 
ceptibility. 



the temperature range 200-300 K, with an effective mo- 
ment jii e ff — 6.1(1) \i b and a paramagnetic Curie tem- 
perature P c± —220(5) K (antiferromagnetic). As will be 
shown in section El the Neel temperature in YbMnOs 
is T/v — 85 K. The large negative P value, such that 
\O p \/Tn — 2.5, could be linked with the frustration of 
the Mn moments on their triangular lattice. Assuming 
for Yb 3+ the free ion value ji e ff = 4.54 we obtain an 
effective moment of 4.1(1) \±b for the Mn 3+ ion (5=2), 
which is lower than the value 4.9 \±b expected for g=2. 
This deviation could be due to our limited experimen- 
tal temperature range where the Curie- Weiss law holds; 
indeed, susceptibility measurements performed in a sin- 
gle crystal up to 350 K 17 obtain the correct effective mo- 
ment with the free ion values for both ions. No anomaly 
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FIG. 3: (Color online) 170 Yb Mossbauer spectra in YbMnOs 
at selected temperatures. The subspectrum corresponding to 
the minority Yb(2a) site (relative occupancy 1/3), is drawn 
in red, that corresponding to the Yb(4b) site in blue. 



in x(T) is found at the magnetic transition, in con- 
trast with magnetic measurements in single crystal s 17 ! 18 , 
where Xc(T) exhibits a small peak at Tn. At low tem- 
perature, a ferromagnetic like increase of x(T) occurs at 
3.5 K (inset of Figj2|), signalling the onset of magnetic 
ordering among the Yb(2a) ions, as confirmed by the 
Mossbauer data to be described next. 

The 170 Yb Mossbauer spectra at selected temperatures 
are shown in Fig.[3j Between 1.35 K and 25 K, they show 
two subspectra with relative weights close to the ratio 
1:2, that we attribute respectively to Yb(2a) and Yb(4b) 
sites. Above 25 K, the spectra are no longer resolved and 
consist of a single line. In Fig. [3l the Yb(2a) subspec- 
trum corresponds to the red line, the Yb(4b) to the blue 
line. At 1.35 K, both subspectra are five-line magnetic 
hyperfine patterns, with a small quadrupolar hyperfine 
interaction. The weakness of the latter precludes any re- 
liable information about the orientation of the hyperfine 
field with respect to the c-axis to be obtained. The hyper- 
fine fields are respectively 117(3) T and 180(3) T for site 
Yb(2a) and Yb(4b). Using the proportionality constant 
C=102T /fi B linking the hyperfine field and the Yb 3+ 
magnetic moment 19 , the saturated moment is 1.15(3) hb 
for Yb(2a) and 1.76(3) fi B for Yb(4b). 

As temperature increases, the Yb(2a) hyperfine field, 
and thus the moment, slowly decreases down to 3 K, 
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FIG. 4: (Color online) Thermal variation of the Yb 3+ (2a) 
relaxation rate extracted from the 170 Yb Mossbauer spectra 
in YbMnOs- The dashed line is a guide for the eye, and the 
red arrow marks the (probably) ferromagnetic transition of 
the Yb(2a) sites. 



reaching 1.00(3) \ib at 3 K. A spectral change occurs be- 
tween 3 and 3.5 K on the Yb(2a) subspectrum: the hy- 
perfine field pattern transforms into a single unresolved 
broad line at 3.5 K (see the red subspectrum at 4.2 K 
in Fig. [3]). Therefore, the hyperfine field vanishes on this 
site above 3.5 K, and the Yb(2a) ion becomes paramag- 
netic. The Mossbauer lineshape then reflects the ionic 
fluctuations within the ground Kramers doublet, and its 
interpretation is postponed until the end of this section. 
The ordering of the Yb(2a) moments below 3.5 K is also 
observed on the magnetic susceptibility, which shows a 
large ferromagnetic-like anomaly at the transition (see in- 
sert of Fig.[2|). The Yb(4b) subspectrum shows no quali- 
tative change from 1.35 K to 25 K: the hyperfine field de- 
creases, reaching 59(5) T (i.e. a moment of 0.58(5) \±b) at 
25 K. At higher temperature, it is not possible to distin- 
guish the Yb(2a) and Yb(4b) sub-spectra. The thermal 
variation of the Yb(4b) moment is presented in section 
IVII in comparison with the moment values derived from 
neutron diffraction. Our 4.2 K spectrum is in agreement 
with that of a recent 57 Fe and 170 Yb Mossbauer study 22 . 

The behaviours of the Yb ions at the two sites are 
thus rather different: above 3.5 K, there is a vanishing 
molecular field at the 2a site, which means that the ex- 
change field from the Mn ions is zero. This is not the 
case for the Yb ions at the 4a sites, which are polarized 
by Mn exchange up to at least 20 K. The neutron study 
described below shows that this exchange field is present 
up to 40 K, and probably up to T/y. This behaviour is 
similar to that in HoMnOs. 

The paramagnetic relaxation of the Yb(2a) ion can be 
interpreted using the Mossbauer relaxational lineshape 
developed in Refl20l. Since the Yb(2a) subspectrum con- 
sists of a single line, the "extreme narrowing" regime 
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holds for the fluctuation rate 7^7, i-e. fe/r ^> A, where 
A is the hyperfine constant associated with the Yb 3+ 
ground Kramers doublet and the excited nuclear state of 
the 170 Yb isotope. It is related to the above mentioned 
proportionality constant C by 21 : A = C gfis gn^n, 
where g and g n are respectively the electronic and nuclear 
gyromagnetic ratios (or g-factors), and ji n the nuclear 
Bohr magneton. Using the Yb(2a) saturated moment 
value 1.15/iB, yielding g=2.3, one gets A — 580 MHz. 
Since the width of the Yb(2a) line narrows as tempera- 
ture increases, the spin relaxation rate increases, accord- 
ing to the expression for the dynamical broadening 20 : 
ATr = 3A 2 r. The fitted values of the relaxation rate -J— 

n Z7TT 

are represented in FigHl It slows down abruptly when 
approaching the transition near 3.5 K from above ("crit- 
ical slowing down"); the T = "extrapolated" value of 
~6 GHz is due to exchange-driven spin-spin relaxation 
and the increase above 10 K can be caused by phonon 
driven relaxation. 
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MAGNETIC STRUCTURE 



The thermal evolution of the NPD patterns is shown 
in Fig. [5j The temperature dependences of the inte- 
grated intensities of typical Bragg peaks (100), (101) 
and (102) are shown in Fig. [6j Strong changes occur 
around Tat = 86 which corresponds to the magnetic 
transition. Below this temperature, an intensity of mag- 
netic origin grows at the (100) and (101) Bragg positions, 
where the nuclear reflections are forbidden by the PG^cm 
symmetry, and the intensity of the (102) peak strongly 
increases due to an additional magnetic contribution. All 
peaks can be indexed within the hexagonal space group 
P6scm and a propagation vector k = 0; this points to 
an AF structure having the same unit cell as the chem- 
ical one, since a ferromagnetic ordering can be excluded 
from the absence of anomaly in the susceptibility at the 
transition. 

The temperature dependence of the (100) and (101) 
Bragg peaks strongly differ from each other. The (100) 
peak remains rather weak down to about 20 K, then 
suddenly increases below. In contrast, the (101) peak 
increases abruptly below 86 K. These variations suggest 
that the (100) and (101) peak intensities are controlled by 
Yb and Mn ordering, respectively, the (102) Bragg peak 
involving contributions of both orderings. Moreover, one 
observes a strong (101) peak in YbMn03, with thermal 
evolution akin to that of the (100) peak in YMn03 16 . 
Recalling that in YMn03, the Mn order does not induce 
a strong (101) magnetic peak, it means that Mn order 
is affected by Y/Yb substitution. Finally, the tempera- 
ture dependence of all magnetic peaks is monotonic, in 
contrast with HoMnOa where Mn magnetic moments re- 
orient at low temperature. 
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FIG. 5: (Color online) Observed and Fullprof calculated NPD 
patterns between 1.5 K and 100 K. The Bragg reflections 
(tics) , and the difference between the observed and calculated 
patterns are plotted at the bottom. 



A. Mn order 

To analyze the magnetic structure we searched for 
all irreductible representations (IR) compatible with the 
crystal symmetry using the theory of group representa- 
tion analysis 23 and the program Basireps 2 -. Description 
of this analysis for hexagonal RMnOs compounds was 
already given in Refflil The atomic positions in the 
unit cell were kept fixed and equal to those determined 
above. In the space groupe P6scm with k = propaga- 
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tion vector, one finds 6 irreducible representations (IR), 
corresponding to the 6 possible configurations of the Mn 
moments reported in Fig. [71 In these configurations, the 
Mn moments lie in the (a,b) plane and their vectorial 
sum is zero. This results from the frustration of the Mn 
moments in their triangular lattice. 

It is important to notice that the I^ and the I^ configu- 
rations are homometric (namely cannot be distinguished 
by NPD) as well as Ti and T 3 25 . The best fit of the data 
was obtained assuming Mn ordering which corresponds 
to the irreducible representation I^ or 1^, in contrast 
with YMnOs where it corresponds to a Ti or 1^ IR. The 
saturated Mn moment value is 3.25 lower than the 
value 4\±b expected for g—2 and S=2. 



B. Yb order 

The magnetic order of Yb moments is more intricate 
and difficult to determine unambiguously using neutron 
diffraction only, since the two crystallographic 2a and 
4b sites order independently. Space group symmetry al- 
lows two types of magnetic orders: along the c-axis or in 
the (a,b) plane, and the coupling between the Yb ions 
can be antiferromagnetic or ferromagnetic (F), leading 
to six possibilities. The Mossbauer spectroscopy data 
show that the low temperature Yb magnetic moments 
are small: 1.76 \±b for Yb(4b) and 1.15 \±b for Yb(2a), as 
compared to the saturated Mn moment value of 3.25 /i^. 
Then the magnetic ordering of the Yb sublattices only 
gives incremental contributions to the NPD spectra, and 
their magnetic structure was solved with the help of 
Mossbauer spectroscopy. 

In the temperature range 3.5 K<T<86 K, we consid- 
ered contributions only from ordered magnetic moments 
on the Yb(4b) sites. For T<3.5 K, we considered also 
ordered moments on the Yb(2a) sites. As to the 4b sites, 
the best agreement with the NPD patterns is obtained 
for Yb moments oriented along the c-axis. We find that 
the Yb moments of the 4b sites are antiferromagnetically 
coupled within a given layer, and that the z and z+^ 
layers are ferromagnetically coupled. At T = 3.7 K, 
this refinement yields a discrepancy factor Rjviag of about 
3.1 %, much better than the values around 20 % given 
by alternative solutions. As shown in Fig.[8j the moment 
values on the Yb(4b) sites deduced from NPD patterns 
are in excellent agreement with those obtained from the 
Mossbauer spectra. 

As to the 2a sites, which order below 3.5 K with much 
smaller moments, we fixed the moments to the values 
deduced from the Mossbauer spectra. The orientation of 
the moments was refined using our NPD datas. There 
are two possibilities for the Yb(2a) to orient : along the 
c-axis or in the (a,b) plane. Taking into account the 
ferromagnetic-like increase of the susceptibility at 3.5 K, 
we assume a F coupling between two adjacent Yb(2a) 
layers. We obtain a discrepancy factor of R=3.05 % for 
the 2a moments in the (a,b) plane, slightly better than 
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FIG. 6: (Color online) Integrated intensities of the (100), 
(101) and (102) reflections between 1.5 and 100K (lines are 
guides to the eyes). 



for the moments along the c-axis (R=3.44 %). It suggests 
that the Yb(2a) moments are perpendicular to the c-axis. 



Finally, the moment values at 1.5 K are 
^Mn=3.23(5)/i J B, ray 6 (4 6 )=1.77(5)/i£ and m Y b{2a) = 
1.13(5) He- The evolution of the moments versus 
temperature is plotted in Fig. [HJ 



6 




Mn spins 



► z=0 
:>z=l/2 



Yb(4b) site 



FIG. 7: (Color online) Symmetry allowed Mn spin orders in 
hexagonal RMn0 3 . In YbMn0 3 , T 2 or T 4 IR yield the best 
fit to experiment. 



VI. DISCUSSION 

In YbMnOa, the Mn magnetic order with or 
I\ symmetry agrees with expectations from optical 
measurements^. Such a type of order is also expected 
from the general tendency observed in hexagonal RMn03 
compounds versus the rare earth ionic radius 26 . T2-4 
magnetic orders seem to be stabilized for low ionic ra- 
dius (Er,Yb,Lu), and Ti_3 for higher ionic radius (Ho, 
Y, Y-Er). This could be related to the amplitude of the 
distortion of the Mn lattice with respect to triangular 
symmetry. 

As concerns the rare earth moments, we found that 
the 46 and 2a sites order independently. This is also the 
case in HoMnOs, where ordered Ho moments on 46 sites 
are observed below 35-40 K, whereas Ho 2a sites remain 
paramagnetic down to 5 K 28 . In HoMnOs, a change of 
the Mn structure (from T2-4 above 38 K to Ti_3 below) 
occurs together with the rare earth ordering. This is not 
the case for YbMnOs, possibly due to the lower value of 
the Yb moment. 

In spite of a very large amount of experimental and 
theoretical work, the nature of the interactions which 
control the magnetic structure in RMn03 compounds is 
still unclea r 27 ! 28 . To shed some light on this point, we 



first compare the measured temperature dependence of 
Mn and Yb moments with calculations from a molecu- 
lar field model. Then we briefly discuss the symmetry 
of interactions necessary to explain the observed mag- 
netic orders. Finally we propose possible origins for the 
magnetic interactions in RMn03, compatible with both 
experiment and calculations. 



A. Molecular field analysis 

The thermal variations of the Mn and Yb magnetic mo- 
ments derived from both neutron and Mossbauer data are 
shown in Fig.[U The Yb(4b) moment values derived from 
both techniques are in remarkable agreement, and their 
thermal variation has an unusual shape, with upwards 
curvature in a large temperature range. We applied the 
molecular field model in order to calculate these thermal 
variations. We first obtained the T dependence of the 
Mn moment by a self consistent mean field calculation 
using the Brillouin function #2 0*0 for S=2 and a molec- 
ular field constant Ao representing Mn-Mn exchange. The 
standard formula has to be slightly modified because the 
saturated Mn moment of m sa t=3.23(5) /i# is lower than 
gSjiB=^ I^b- So we used the following expression: 



m Mn (T) m sat B 2 



gjjL B S Xom M n(T) 



(1) 



which ensures that saturation occurs at m sa t=3.3 \±b- 
We obtain a good fit of the experimental data (blue 
line in Fig. [8]) with A =19 T/fi B . Then we fitted the 
T dependence of the Yb(4b) moment assuming that the 
ground crystal field doublet alone is populated, i.e. de- 
scribing the Yb ion by an effective spin 5=1/2, with a 
g-factor gyb m the direction of the net exchange field 
from the Mn ions. Since the saturated Yb(4b) moment 
is mo = \gyb^B = 1-75 he, then: g y6 =3.5. The Yb(4b) 
ion is submitted to a molecular field XirriMn from Yb-Mn 
exchange, since we neglect Yb-Yb exchange. Its moment 
is then calculated using the expression: 



m Y b(T) -gvbl^B tanh 



gvbl^B \im M n(T) 



2k B T 



, (2) 



where rriMn(T) is obtained from expression (pQ). We ob- 
tain a good fit to the data (magenta line in Fig. [8]) with Ai 
= 2.1 T/fiB- The unusual upwards curvature comes from 
the tanh function with a small argument, the constant Ai 
being much smaller than Ao- With this model, a small 
ordered Yb moment should persist up to Tat = 86K, with 
a value below 0.3 beyond the accuracy of neutron 
and Mossbauer probes. We conclude that the Yb ions 
on the 46 sites orient in a non-vanishing net molecular 
field through Yb-Mn exchange. As to the Yb moment 
on the 2a sites, their much lower ordering temperature 
(3.5 K) strongly suggests that they orient in their own 
molecular field, through Yb-Yb interactions. From this 
analysis follows the hierarchy of exchange interactions in 
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FIG. 8: (Color online) Observed Mn, Yb(4b) and Yb(2a) mo- 
ments from NPD and Mossbauer measurements (dots) versus 
temperature. The solid lines correspond to the molecular field 
calculation. 



YbMn03, in decreasing order : Mn-Mn, Yb-Mn and Yb- 
Yb. Similar conclusions could be valid for other R ions 
(Ho, Er) in hexagonal RMnOs. 



IR 


A-pair 


R(4b) coupling 
Intraplane Interplane 


Ti 


(+ -) or (- +) 


AF AF 


r 2 


(" ") 


F F 


r 3 


(- +) or (+ -) 


F AF 


r 4 


(+ +) or (- -) 


AF F 



TABLE II: R(4b) magnetic orders deduced for each configu- 
ration Ti of the Mn moments: i) by our analysis of the local 
magnetic environment using the sign of the A-pairs; ii) by 
the minimization of DM and dipolar energies. Both methods 
yield the same result. 



in environments with identical (resp. opposite) A-pairs 
align parallel (resp. antiparallel) along c, one can pre- 
dict the R magnetic order on a 4b site for a given Mn 
order. This is reported in Table HTl This description pro- 
vides a simple way to distinguish two homometric con- 
figurations, knowing the order on the 4b sites. Experi- 
mentally, the intraplane Yb(4b) coupling is AF while the 
interplane coupling is F. Therefore, the Mn order must 
be of T4 type and T2 can be discarded. This description 
also agrees with recent experiments^*^! in HoMnOs and 
can be easily extended to the rest of the RMn03 family. 



B. Mn environment of a R ion 

Assuming that R(4b) moments order via R-Mn inter- 
actions, we consider the Mn environment around a R(4b) 
site, taking into account nearest neighbor Mn ions only. 
In the following we will consider only the first four IR 
(Ti to T^) which describe the magnetism of the whole 
RMnOs family (see fig 5). The moment on a R(4b) site 
is expected to orient along the c-axis, to obey the trigonal 
symmetry of its Mn environment. This agrees with our 
observations in YbMnOa and with the magnetic order in 
HoMnO^^. 

In order to depict the local magnetic environment of 
each R(4b) site, we need to take into account the orienta- 
tion and position of the Mn moments. For this purpose, 
we define the A factors as : 

Ai = ^^Sj.fij^ (3) 

3 

where the sum runs over the three Mn first neighbors of 
the R(4b) site in a particular plane, Sj is the magnetic 
moment of Mn at site j and fij is the vector joining Mn 
site j to R(4b) site i. 

The A factors are either positive (+) or negative (— ) 
for a given Mn plane depending of the magnetic config- 
uration T^ of the Mn ions. This leads to four possible 

magnetic environments for the R(4b) ion, (+ +), ( ), 

(+ -) and (- +), defined by the signs of the (A" p Af° wn ) 
pairs associated with upper and lower Mn planes. One 
can easily see that there is either only one or two op- 
posite A-pairs allowed for each IR, whatever the IR re- 
tained for the Mn order. Considering that R moments 



C. Origin of the molecular field 

The exchange energy between a Yb ion and its Mn 
first neighbors is classically written as a sum over Mn 
sites JSMn-Sybi where J is the isotropic first neighbor 
exchange interaction. Each Yb atom (either on a 2a or 
4b site) is at the top of a pyramid with an equilateral 
triangular basis of Mn ions. Therefore, the exchange field 
induced at either Yb site cancels by symmetry for any 
of the Ti_4 representations of the Mn magnetic order. 
Moreover, the classical exchange term cannot explain the 
perpendicular orientations of the R and Mn moments. 

Experimentally, the exchange field is effectively zero at 
the Yb(2a) site (above 3.5 K), but it is not at the Yb(4b) 
site. So, in order to explain the R(4b) magnetic order, we 
have to introduce an interaction which is not isotropic, 
i.e. which depends on r^-, the radius vector linking a 
R(4b) to a Mn ion. This is the only way to obtain an out 
of plane effective field. 

In the following we consider two interactions fulfill- 
ing this condition: the magnetic dipolar interaction and 
the Dzyaloshinskii-Moriya (DM) coupling between R(4b) 
and Mn spins. The magnetic dipolar interaction is de- 
fined as : 

7~l>dip = —^^(JiBij, (4) 

where Bij is the dipolar field induced by Mn moment j 
at R(4b) site z, and Si the spin of the R(4b) ion at site i. 
The Dzyaloshinskii-Moriya (DM) interaction between 
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the Yb(4b) and the Mn moments is written as 

n DM = ^Dij.fa x Sj). (5) 



According to the crystal symmetry and the rules defined 
by Moriya 32 , Dij is perpendicular to the mirror plane 
including both R and Mn sites, i.e. to the c-axis, and 
defined as = r]cx , where r] is a spin-orbit coupling 
constant. 

Although the dipolar and DM interactions are intrin- 
sically different, they may be expressed in similar ways 
in terms of the A factors: 

H dip = J2^i( A i rZ T P + ^ r z ) down ] (6) 

h dm = nY^vMT + A d ™ n ). (?) 



Minimizing the energy for these interactions should 
therefore lead to the same configurations as found in sec- 
tion B. To confirm this, we have rewritten both interac- 
tions in terms of a matrix J^— such as TL = ■ Si Jij Sj , 
where is expressed as: 



Ji 



DM 





—D z 

1-3 

D y . 
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Df 



jdip. 



—D x 

•J 

ij ij 
'ij'ij /;; 



ij ij 
ij ij 

v ij ij ij ij ij ij / 




(8) 
(9) 



It is worth noting that the above matrixes both have 
non-diagonal terms, which depend on the R-Mn pair con- 
sidered. This is a necessary condition to ensure a non- 
zero interaction between Mn moments in (a,b) planes and 
R(4b) moments along the c-axis. 

We determined the stable Yb(4b) configuration for 
each IR of the Mn order. We find that both dipolar and 
DM interactions give the same result, in agreement with 
our analysis of the magnetic environment made above 



(see Table HI]) . Namely they stabilize the T4 configura- 
tion of the Mn moments for the Yb(4b) order determined 
experimentally. 

The analysis described here provides the conditions re- 
quired for the magnetic interactions to yield the R and 
Mn magnetic orders observed experimentally. To go fur- 
ther, a quantitative calculation of the energy should be 
made. Our estimations made for the the dipolar inter- 
action show that it is too weak by one order of mag- 
nitude to fully account for the molecular field constant 
^Mn-Yb(=^i) = 2.1 T\i^ found in section A, Fig. [HJ As 
for the DM interaction, an estimation would require an 
evaluation of the relevant spin-orbit coupling constant. 

We hope that this analysis will stimulate further 
theoretical work on the R-Mn interactions at the origin 
of the molecular field. This will allow one to understand 
the microscopic origin of the rare earth magnetic order 
in the hexagonal RMn03 family. 

VII. CONCLUSION 

We have studied the magnetic order in YbMnOa and 
found a consistent description by combining Mossbauer 
and neutrons probes. The temperature dependences of 
the Yb and Mn magnetic moments are well fitted in a 
molecular field approach, showing that Yb(4b) moments 
order due to the Mn molecular field, whereas the Yb(2a) 
moments order at much lower temperature through Yb- 
Yb interactions. 

We could determine the Mn magnetic structure and 
distinguish between two homometric configurations by 
considerations about the magnetic environment of each 
R(4b) sites. In YbMn03, we found a I\ configuration for 
Mn moments associated with F stacking of AF Yb(4b) 
layers. This approach can be extended to the whole 
hexagonal RMn03 family. We propose possible mecha- 
nisms for the R-Mn interactions to be studied by further 
theories. 
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